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Abstract 

We derive conservation and balance laws for the translational gauge theory of 
dislocations by applying Noether's theorem. We present an improved translational 
gauge theory of dislocations including the dislocation density tensor and the disloca- 
tion current tensor. The invariance of the variational principle under the continuous 
group of transformations is studied. Through Lie's-infinitesimal invariance criterion 
we obtain conserved translational and rotational currents for the total Lagrangian 
made up of an elastic and dislocation part. We calculate the broken scaling current. 
Looking only on one part of the whole system, the conservation laws are changed 
into balance laws. Because of the lack of translational, rotational and dilatation 
invariance for each part, a configurational force, moment and power appears. The 
corresponding J, L and M integrals are obtained. Only isotropic and homogeneous 
materials are considered and we restrict ourselves to a linear theory. We choose 
constitutive laws for the most general linear form of material isotropy. Also we 
give the conservation and balance laws corresponding to the gauge symmetry and 
the addition of solutions. From the addition of solutions we derive a reciprocity 
theorem for the gauge theory of dislocations. Also, we derive the conservation laws 
for stress-free states of dislocations. 

Keywords: Conservation laws; balance laws; Noether symmetries; dislocations; 
gauge theory. 
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1 Introduction 



In the Lagrangian field theory of elasticity the existence of inhomogeneities in an elastic 
medium is expressed through the explicit dependence of the Lagrangian on the position. 
This shows up as a configurational force acting on the inhomogeneity being an elastic 
singularity. Understanding elasticity as a classical field theory, Eshelby [H, 0| and Morse 
and Feshbach 0] derived the elastic energy-momentum tensor, called also Eshelby stress 
tensor. Supplementary, Eshelby [lj and Morse and Feshbach 0] introduced independently 
the field intensity vector and the field momentum vector of elasticity. 

Later Giinther 0] used the Noether theorem to find all the continuous transformation 
groups leaving the potential energy density of linear elastostatics invariant. With the 
known Lie-point symmetries he constructed conserved currents and path-independent 
integrals and related them to torsion and bending problems of rods, plates and shells. 
Knowles and Sternberg Q extended these results for the case of finite elasticity. Fletcher 
[sl continued further and calculated the variational symmetries and conserved integrals in 
linear elastodynamics. Moreover, Olver |7| investigated the Lie-point and Lie-Backlund 
symmetries of elastostatics and derived the corresponding conservation laws. An overview 
of conservation and balance laws of elasticity can be found in the books of Maugin and 
Kienzler and Herrmann j^]. 

In the gauge theory of classical electrodynamics the inhomogeneous Maxwell equa- 
tions can be obtained from a Lagrangian density after varying the scalar and vector 



gauge potentials. The Lagrangian density consists of two parts [10|, [III], describing the 
electromagnetic field and another one being responsible for the interaction between the 
electromagnetic field and the charge and current density. Kadic and Edelen 12| proposed 
a theory of dislocations in a continuum with gauge invariant equations of motion. The 
formal structure of these equations is analoque to the inhomogeneous Maxwell-equations. 
Kadic and Edelen [l| and Edelen and Lagoudas [HJ presented a mathematical gauge 
theory of Yang-Mills type for elasticity including dislocations. They derived equations of 
motions and suggested some improvement since the given constitutive laws were inade- 
quately chosen. In general the formal structure of the theory was correct, since canonical 
energy-momentum tensors and configurational forces such as the Peach-Koehler force 15 
has been successfully calculated. 

In this paper we give a complete system of equations of motion for the translational 
gauge theory of dislocations. We calculate the Lie-point symmetries of this system and 
also the variational and divergence symmetries of the Lagrangian, following the standard 



technique of Lie-group analysis given in the book of Olver [16|j. The canonical energy- 
momentum and angular-momentum tensors, as well as the dilatation current vector are 
determined. The generalization of the conserved integrals from the linear theory of elas- 
todynamics to a translational gauge theory of dislocations is given. Symmetry-breaking 
follows by taking into account only the elastic or dislocation part of the system. In this 
case we calculate the configurational force, moment and power acting on dislocations. 
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2 Gauge theory of elasticity 



According to the field theory of elasticity, the equations of motion for a compatible con- 
tinuum can be derived from a variational principle with the Lagrangian of linear elasticity 

1 1 

C e = T e - W e = -piiti - - Oij Uij. (2.1) 

The canonical conjugate quantities to the kinematical velocity iii and displacement gra- 
dient U{ j £11*6 thereby defined as follows 

Pi '= W '= (2 ' 2) 

These are the momentum vector and the force stress tensor a^. Since the works of 



Kadic and Edelen [12j, [13| and Edelen and Lagoudas [14| it is known that symmetry 
breaking of the homogeneity of the action of the three-dimensional translational group 
T(3) in elasticity leads to a gauge theory of dislocations. In terms of symmetry, we require 
from the Lagrangian density (12.11) to stay invariant under the internal transformation of 
the displacement field 

u!i = Ui- ft (2.3) 

given with an arbitrary vector-function fi(t,x). The transformation (12.31) represents the 
generalization of a rigid body translation with = const. This requirement can only be 
fulfilled by introducing the gauge potentials (p and 4>. They restore the invariance of (12.11) 
through the following two steps. First of all, a minimal replacement 

Pij ■= ^ jUi = Uij + (pij , Vi := V t Ui = Ui + (fi, (2.4) 

of the displacement gradient u^j and Newtonian kinematic velocity Ui by the translational 
gauge-covariant derivatives VjUi and V t Ui is used. These represent the elastic distortion 
tensor and elastic velocity vector Vi, respectively. Both quantities are now incompati- 
ble. The gauge potentials may be identified with the plastic velocity and plastic distortion: 
ipi = —vf and 4>ij = —0ij. On the other hand, for compensating the time fi and space 
fij derivatives in the Lagrangian (12. ip caused by (12.31) the gauge fields have to transform 
as follows 

¥'i = Pi + fi, o{j o,j ■ f,. r (2.5) 

The transformations (12. 5p and (12.31) are making up the local gauge transformations leaving 
the Lagrangian density 

1 1 

C c = -p^- -a^P.j, (2.6) 

form-invariant called gauge-invariant. The elastic Lagrangian density (12. 6p describes an 
incompatible elastic material. By construction a local gauge symmetry is valid 

v' t = v u [3^ = /%. (2.7) 
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Dislocations being now present contribute also by themselves to the total energy of the 
system. With the gauge fields ip and <p we can introduce the two translational field 
strengths 



Tijk 



I i i 



or in terms of (3^ and Vi 



Tijk Pik,j l^ij,ki 



-(fij + (f>ij 



(2i 



(2.9) 



called the dislocation tensor (torsion tensor) and the dislocation current tensor, respec- 
tively. These are the kinematical quantities of dislocations. They are related to each 
other through the Bianchi identities 



£jkiTijk,i — 0, T^k + 2 1{[j,k] — 0. 



(2.10) 



The transformation of the fields (iij, ip^ 0^) according to the local gauge transforma- 
tions (12.51) and (12.31) leaves the quantities 



T' 



Ta 



ijk i 



gauge-invariant. The torsion tensor was introduced by Elie Cartan 



eralization from Riemannian to non-Riemannian geometry. Kondo |18|, [I9j, Bilby et al 
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(2.11) 



as a gen- 



[20j, |21] and Kroner 22] understood the meaning of the three-dimensional torsion tensor 
as the dislocation density tensor. The Lagrangian of the dislocation field can be written 
in the following form 



C 



di 



T rii - W, 



di 



2 T)ijlij 



^ T^ijkTijki 



(2.12) 



where and W& describe the kinetic and potential energy density of dislocations. The 
total Lagrangian density describing the whole system (in linear approximation) is given 
as the sum of (j2SD and (l2~T2l) 



C = T-W = C e + C 



di 



^PiVi 



&ijftij ~t~ T^ijlij 



The canonical conjugate quantities can be obtained from C as 
dC ^dC „ dC 



Pi ■-- 



i.i 



dhj ' 



^ H-ijkTijk- 



dC 



m. 



tjk 



[2.13) 



(2.14) 



Pi, a^, Iij, and are the momentum vector, the force stress tensor, the dislocation 
momentum flux tensor, and the pseudomoment stress tensor, respectively. 

The Euler-Lagrange equations derived from the total Lagrangian C = C{vi, fy, Zy, r^ fe ) 
are given by 



D dC p dC 

* diii 3 duij 

n dC ^ dC 
D t — + Dj 



dpi 



Fir 

EfAC) = D t -^f - + D fe 



dpij 
dC 



0. 

dC 

dpi 
dC 



'10 



(2.15) 
(2.16) 
(2.17) 
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where T) t and Dj are thereby the so-called total derivatives: 



d <9 d 

V t = 1 ^+u al — + <p al — + <j) a p 1 r— + ... (2.18) 

Dl = ^ + Ua * ^ + Ua " d^~- + ^ W + 0a/M W~« + • • ' • (2 ' 19) 

Written in terms of the canonical conjugate quantities (12.141) . Eqs. (I2.15I) - (I2.17I) take the 
form 

D tPl - DytTy = 0, (2.20) 

D j D ij +p i = Q, (2.21) 

D t D y + D fc i% + o-y = 0. (2.22) 

Eq. ( 12.201) represents the balance of linear momentum. Eqs. ( 12.211) and ( 12.221) are the 
equations for the balance of dislocations. The coupled system (I2.20j) ( T2.22j) together with 
the Bianchi identities ( I2.10p are analogous to the inhomogeneous Maxwell equations. Like 
in electromagnetic field theory by coupling the field to external charges and currents (p, j), 
the dislocation field will be influenced by the momentum and stresses (p, cr). The force 
stress cr and linear momentum p works like sources driving dislocations. The elastic 
medium plays also the role of transmitting the interaction like an elastic aether. The 
conservation of linear momentum appears as an integrability condition for the balance of 
dislocation equations. This can be seen by applying D t on (12.211) and on (12.221) and 
subtracting the last from the first one. It plays the same role like the conservation for the 
charge in the electromagnetic field theory. The most general linear constitutive equations 
for the momentum, the asymmetric force stress, the dislocation momentum flux tensor 
and the pseudomoment stress of an isotropic and centrosymmetric medium are 

Pi = pvu (2.23) 

Oij = \5ijPkk + vifaj + (3ji) + -yifiij - pji), (2.24) 

Dij = dtSijlkk + d 2 (Iij + Iji) + d 3 (I i:j - Iji), (2.25) 

Hijk = c{Tij k + c 2 (Tj ki + T ki j) + c 3 (8ijTu k + S ik Tiji), (2.26) 

where p is the mass density and with 9 material constants fj,, v, 7, ex, . . . , c 3 and d±, . . . , d 3 . 
Here we have used a constitutive relation for the dislocation momentum flux tensor 
Da = L 

ijkdki with Lij k i — L k uy A similar constitutive law was used for material with 



microstructure by Maugin [23j. Maugin [23] also noted that for the isotropic microinertia 



is often assumed a too simple form which is not imposed by the formulation of the theory. 
In Eq. (12.251) we have assumed the most general one for an isotropic theory. On the other 



hand, Kadic and Edelen [13J and Edelen and Lagoudas |l4j used the constitutive relation 
= S2 Iij which is to simple. 

The requirement of non-negativity of the energy (material stability) E = T + W > 
leads to the conditions of positive semidefiniteness of the material constants. Since Vi, 
Pij, 1^ and Tij k are uncoupled from each other, the conditions can be studied separately: 
T e > 0, T(ji > 0, W e > and W^i > 0. Thus, the characteristic constants of the material 
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have to satisfy the following conditions 

P>0, 

d 2 > 0, d 3 > 0, 3dx + 2d 2 > 0, 

/i>0, 7>0, 3A + 2/i>0, 

ci - c 2 > 0, ci + 2c 2 > 0, ci - c 2 + 2c 3 > 0. (2.27) 

Substituting the constitutive equations in the above system (I2.20I) - (I2.22I) . we get for 
the field variables Ui,(j>ij,ipi a system of 15 coupled linear partial differential equations 
A = (A 1 ,...,A 1B ) = 0: 

p(ui + ipi) - \{u jtji + <f> jjti ) - (/j, + "i){u idj + (pijj) - (fj, - i)(u jtij + <pji,j) = 0, (2.28) 
diifjji ~ 4>jj,i) + (d 2 + d 3 ){<p idj - 4>ij d ) + {d 2 - d 3 )(tp jAj - (fijij) - p(ui + ipi) = 0, 

(2.29) 

di5ij{(p k ,k - (f>kk) + (d 2 + d 3 )((p id - (fiij) + (d 2 - d 3 )(ip jti - 0^) (2.30) 

— Cl(<j)ik,jk — (fiij,kk) ~ c l{$ji.kk — 4>jk,ik + <Pkj,ik ~ ^kijk) ~ C 3 [8ij(<filk,lk ~ <fiu,kk) 

+ {4>kk,ji - 4>kj,ki)] ~ ^iji u k,k + 4>kk) - + i){ui,j + 4>ij) -((*- i){v>j,i + 4>ji) = o. 

Eq. (12.281) is a generalized inhomogeneous Navier equation for u. Eq. (12.291) has the form 
of a generalized inhomogeneous Helmholtz equation for ip and Eq. (I2.30p is a kind of gen- 
eralization of an inhomogeneous Klein-Gordon equation for <f>. Due to the inhomogeneous 
parts they are coupled. 

3 Lie symmetries 

The infinitesimal continuous transformation acting on the independent (t, x) and depen- 
dent (u, (p, 4>) variables build a Lie group G. If G is the Lie group of invariance of the 
system (I2.20I) - (I2.22I) . then it is also a symmetry group. The infinitesimal group action for 
the independent and dependent variables has the form 



x- = Xi + eXi(x,t,u, (p,(p) H (3.1) 

t' = t + er(x,t,u, tp, (f>) H (3.2) 

u' a = u a + EU a (x,t,u,(p,<j>) H (3.3) 

fa = + £^ a (x,t,U,ip,(f>) H (3.4) 

0a/3 = <Pal3 + £ Qaflfa, t,U, (fi, (/))+■■■ , (3.5) 
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where the infinitesimal generators are defined by 

Xi(x t t,u,(p,<f>) := 
r(x t t,u,(p,<f>) := 
U a (x,t,u,(p,(f>) := 

<5> a p{x,t,u,Lp,c/)) = 



dx'j 

de 

dt' 

de~ 

du' 



e=0 



E=0 



de 

d<f' a 



de 



e=0 



e=0 



de 



(3.6) 
(3.7) 
(3.8) 
(3.9) 
(3.10) 



e=0 



These infinitesimal generators build the following vector field 

d_ 

dt 



9 v d TT d d _ d 

V = T — + Xi — + U a — + ^ a + Tr — 



dxi 



du a 



dip c 



(3.11) 



>afi 



Since the system (I2.28l) -( |2.30l) is of second order, we need the prolonged vector field of 
second order 



pr^ 2 V = v + a + b + c, 



where the vector fields a, b, c read 
a = U, 



du, 



9 FT 9 FT 

+ u at — + u ( 



b = ty 
c — $ 



d 



du a 
d 



a/j 



d - d - d 

+ Uait TT. \~ Uatt 



^ at oi • oiij 



du a ij 
d 



du a . 



du r 



dipa 



dip, 



d(fa,i 

d - d - d 

777 ^ "z 3 * 771 ^ 



-d-d 

+ ^ait 77— + * att 



dip 



■>a/3,i 



d(f) a 



®aj3it 



dip. 

d 



d 



a/3tt 



(3.12) 

(3.13) 
(3.14) 
(3.15) 



with 



(3.16) 
(3.17) 
(3.18) 
(3.19) 
(3.20) 

Similar expressions are given for the coefficients of the vectors b and c if in the above 
formulas the variables U a and u a are substituted by ip a and $ a/ g, (p a p, respectively. 
The Lie-group G is a group of invariance of the system A if and only if |l6j 





= D;(*7 a - 


XkUa,k — 


TU a ) + X k U aM + TUaj, 


Uat 


= V t (U a - 


X k U a ,k — 


TU a ) + X k U a ^ k + TU a , 


U a ij 


= DJVK* 


— X k u ak 


T^a) X k U a k ij ~\- TU a {j, 


U a it 


= DiD t (U a 


— XkU a ,k 


- TU a ) + X k U ajk i + T<il aji , 


U a tt 


= VMU« 


— X k u ayk 


- ru a ) + X k u a>k + r u a . 



pr^ 2 V(A) = 0, whenever A= 0, 



(3.21) 
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for every infinitesimal generator v. Applying the procedure (13.211) to the system of 
Eqs. (I2.28p -( l2.30p . we get for the infinitesimal generator i> 



p{U m + ty it ) - Ml',,, + <$> j:ji ) - (// • 7)(%? + - !/' ' iWjji + 



A=0 



A=0 



(3.22) 

= 0, 
(3.23) 

diSiji^kkt - $ kktt ) + (d 2 + da)(Vijt ~ hjtt) + {d 2 - d 3 )(^ jit - $ m ) 

~ Cl(^ikjk ~ ^ijkk) — c 2^jikk ~ §jkik + §kjik ~ &kijk) ~ C3 {&lklk ~ &llkk) + i^Uji — 



\5ijiUu + $„) - (/i + 7 )(E/i i + $y) - - t)(^ + %») 



0. 



A=0 



From the condition above the form of the infinitesimal generator i> is found 

r = d 4 , 

Xi di ~\~ 6ij k (ljX k} 

U a CU a ~\- ^aj k OijU k fa ^ai 

ty a = ctp a + e ajk aj(p k + f a + (p a , 
$ Q/ g = c0 a/ g + e Q j k aj(p k p + epj k aj<j) ak + f a ,/3 + 4>ai3, 



(3.24) 



(3.25) 
(3.26) 
(3.27) 
(3.28) 
(3.29) 



where a,, di and c are arbitrary parameters. Moreover, f a is a gauge function and u a , tp a , 
(pap are solutions of the Euler-Lagrange equations (I2.20I) - (I2.22I) . The symmetry algebra 
is generated by the following vector fields: 



Of 

A 

dxi 



d 

9 a. 
v =Ui- h (Pi 

OUi 



+u ^ k + ^ 



d d d 



?/,:/ 



(translation in time) 

(3.30) 

(translations in space) 
(3.31) 

(rotations) (3.32) 



9 j. 

d(pi 



d 



5 - f ^ 4- f ^ f 9 
V Ti r Ji ~7{ r Ji j 7T7 

OUi Ofi a( Pij 

6- a -9 — d 
v = Ui — h (Pi ttt + <pij ^— 



(scaling) 
(gauge) 



(3.33) 
(3.34) 



dui 



dipt 



(addition of solutions). 

(3.35) 
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4 Conservation and balance laws 



Since the famous theorem of Noether [24| , it is well-know that to each of the continuous 
symmetries of the Lagrangian a conservation law for the physical system corresponds. 

The Lagrangian C depends on the first derivatives of the dependent field variables. 
In that case the infinitesimal criterion of invariance [l6.] says that a Lie group G is a 
variational or divergence symmetry of C if and only if 

pr«i/(£) + C (D t X t + D t r) = D t B t + D t B A , (4.1) 

where Bi and B 4 are opportune analytic functions. If Bi ^ and B4 ^ 0, then i> is 
the generator of a divergence symmetry of the Lagrangian C. If B 4 = and Bi = 0, u 
generates a variational symmetry of C Every variational or divergence symmetry of the 
Lagrangian C is also a symmetry of the associated Euler-Lagrange equations. The inverse 
statement is not true. 



After some standard calculations (see, e.g., Lazar [251 ] ) we find the conservation law 
in characteristic form 

+ D t A, + QlEl{C) + QK(£) + Q%E%{C) = 0, (4.2) 
where the characteristics are defined by 

Ql = U a - XjU aJ - TU a , (4.3) 
Ql = ^! a - Xjtpaj - T <p a , (4.4) 

Qtp = $ "/3 ~ Xj'fiaPJ ~ T< fial3- (4.5) 

Therefore, if the Euler-Lagrange equations fl2.15p - fl2.17p are fulfilled, we speak of a con- 
servation law 

D t A A + DiAi = 0, (4.6) 

where Ai is the associated flux and A 4 is the conserved density. By the divergence theorem 
one finds the conservation law in integral form 



D t A 4 dV+ / Ain t dS = 0. (4.7) 
v Js 

In case where the equality (14.61) is not fulfilled we speak of a balance law 

D t A A + DiAi ^ 0. (4.8) 

4.1 Canonical currents 

If a variational or divergence symmetry is given, then the corresponding components A 4 
and Ai of the conservation law (14. 6 j) are given by 

bc bc bc 
* = £T + <B _ + g S _ + Qj, ijs --*, (4.9) 

BC BC BC 
A l = CX l + Ql — + Ql — + Q* — Bi. (4.10) 
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4.1.1 Translational invariance 

The translation acts only on the independent variables. The Lie-point group transforma- 
tion of the translation in space and time is given by the formulas 

x[ = Xi + e k 5 ki , t' = t + e 4 5 44 , (4.11) 

leaving the field variables unchanged 

u' a = u ai (p' a = ip a , 4>' af5 = (j) af3 . (4.12) 

The components of the generator (13.111) corresponding to the infinitesimal transforma- 
tions (14.111) and (14.121) take the form 

X ki = 8 ki , t = S u , U a = 0, V a = 0, $ a/3 = 0. (4.13) 

Using Eqs. (14.91) . (14.101) and (14.131) . we obtain for the translational density and flux 
quantities 

a r r ®£ dC ± dC , 

A ki = Ld ki - U a , k - (p a k 0a/3,fc7T7 , ( 4 -14) 

du a>i dip a>i 0(p af 3,i 

A k4 = -U ak — p a ,kTT-- <Pal3,k—; , ( 4 -15) 

OU a 8if a d(p al3 

dC dC ■ dC . 

A u = -u a - (p a - Papwz , (4.16) 

OUa,i 0<fa,i 0<Pa/3,i 

. 9jC . dC ■ d£ 

A 4i = L-u a — f a — Pap—- — • (4-1 7 ) 

du a dip a dcf) a p 

In terms of the momentum p^, the dislocation momentum flux Dij, the force stress err- 
and pseudomoment stress they read as 



Pi 



ki 



v k 

Si 
H 



—A ki — —C 5 ki — <J a i U a ,k — Dai Pa,k + Hafii 4>a/3,k, (4-18) 

Ak4 = ~Pa Ua,k - D a/3 (f) a p,k, (4.19) 

A 4 i = Oai Ua + D ai tp a - H af}i <j) a p, (4.20) 

-A A4 = -£ + p a u a + D a/3 4> aj3 . (4.21) 



The tensor P k i is the canonical Eshelby stress tensor of dislocation gauge theory. The 
vector V k is the canonical pseudomomentum or field momentum density and the vector 
Si is called the canonical field intensity or Poynting vector. The scalar Ti is the canonical 
energy density. They are related to each other by the following local conservation laws 

D t V k - = 0, (4.22) 

D t H - DA = 0. (4.23) 

The first equation represents the canonical conservation law of pseudomomentum, while 
the second one constitutes the canonical conservation law of energy for the dislocation 
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gauge theory. By the help of the Gauss theorem, these conservation laws appear in integral 
form 



Jk 
I 



P kini dS - / D t V k dV = 0, 
s Jv 

:= [ SiHidS- [ D t HdV = 0. 
Js Jv 



(4.24) 
(4.25) 



4.1.2 Rotational invariance 

The three-dimensional rotations act on both, the independent and dependent variables. 
The infinitesimal transformations of the Lie group action are given by 



X^ X{ "T ^ijk^- j^kt 
t' = t, 



U, 



<£' a = <-Pa + tajkZj^k, 
4>'af3 = $oifi + (£ajk4>kf3 + £f3jk(fiak)£j > 

from which the corresponding components of the generator (13. lip are calculated 



X 



it 



Substituting these components into the Eqs. (14.91) and f)4.10p . we obtain 



in i 



dC 



ou. 



dC 

d(p a . 



dC 



dC dC , dC dC 

+ U m h (p m h (Pml — V <Plm 



du jti d(p jti 



Ak4 — £k 



dC 



dC 



dC 



dC 



a. 



du 



— + ¥m TTT- + (pml —■ h <Plm — f 
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dC dC 1 dC 



d(p c 



a/3 



(4.26) 
(4.27) 
(4.28) 
(4.29) 
(4.30) 



£ijk%k, U a j — € a j k U k , ^> a j — tajk^Pk, &a/3j ~ ^ajk4>kp + tpjk&ak- (4-31) 



(4.32) 



(4.33) 



We can introduce the total canonical angular-momentum tensor and the material inertia 
vector according 

M ki : = — A ki = ekmj (x m Pji + <7ji u m + Dji (p m — Hju (p m i — Hiji (f)i m ), (4.34) 
M k : = A M = e kmj (x m Vj + Pj u m + D jt 4> ml + D l} 4> lm ). (4.35) 

It is possible to decompose the total canonical angular-momentum tensor and the material 
inertia vector into two parts 



M kl = M<? } + M, 



0) 
ki ' 



M k = M ( k o) +M k l \ 



(4.36) 
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where 

M u =e kmj x m P ji , (4.37) 

M$ = e kmj (o-ji u m + Dji (p m - H jU cj) m i - Hiji lm ), (4.38) 

M ( ° ] = e kmj x m V 3 , (4.39) 

Mj? = e kmj (pj u m + <j) ml + Dij (f> lm ), (4.40) 

are the orbital and intrinsic (or spin) angular-momentum tensors and material inertia 
vectors, respectively. For the balance of the material inertia vector we obtain 

D t M k - D t M kl = D t M { k o) - DiM® + D t M$ - D,M«. (4.41) 

Due to the conservation of the pseudomomentum (I4.22p . the orbital parts become 

D t M k ° ] - DiM® = e kjm P jm . (4.42) 

Because of the rotational invariance, by using the minimal replacement (12.41) and the 
Euler-Lagrange equations (I2.20l) - fl2.22l) the balance law of the intrinsic angular-momentum 
becomes 

VtM® - DjM^ = e k mj{vmPj - crjiPmi + Djil ml - ]- H ja T mil - H Uj T Um (4.43) 

With Eqs. (I4.18p . (I4.42p and (I4.43P the rotational balance of the total angular momentum 
reads 

Dj-Mfc — DjM fc j = € krn j(v m Pj — Pml^jl ~ Plm^lj + ImlDjl + hmD\j ~ ^ TmilHjil ~ TiimHiij). 

(4.44) 

Using the relations for the dislocation density tensor T ik \ = e k ijOiij and pseudomoment 
stress tensor H ik i = e k ijHij, we obtain 

Di-Mfc — DjMfcj = £kmj(v m Pj — Pmldji — film&lj + ImlDjl + hmDlj — a m [Hji — a\ m Rij) . 

(4.45) 

In the case when the material is isotropic, by using the constitutive Eqs. (I2.23l) -( l2.26l) . 
it can be shown that f)4.44p becomes a conservation law. Using the Gauss theorem, we 
introduce the dynamical X-integral of dislocation gauge theory: 



L k := / M ki mdS- / B t M k dV. (4.46) 
<s Jv 



4.1.3 Scaling 



In the gauge field theory of dislocations the scaling group (dilatational group) is not a 
variational symmetry. The self-similarity is broken, since the defects like dislocations 
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in the material introduce characteristic length scales. Still it is possible to calculate the 
broken scaling quantities. The scaling group acts in infinitesimal form on the independent 
and dependent variables in the following manner 



/ _ 


+ x h 


(4.47) 


t' = 


(1 + e) t, 


(4.48) 


u' a = 


(1 + ed u ) u a , 


(4.49) 


V>'a = 


(l + ed^tpa, 


(4.50) 


<JaP = 


(1 + ed<p) fap. 


(4.51) 



These transformation determines also the form for the components of the infinitesimal 
generator 

Xi = Xi, t = t, U a = d u u a , ^ a = dtpifa, $ Q/3 = d^fap, (4.52) 

where du^d^^d^ denote the (canonical) dimensions of the displacement vector u and the 
gauge potentials ip, </>. If one substitutes the relation (14.521) into the Eqs. (14. 9[) and (14.101) . 
one obtains for the scaling quantities 

dC dC 

Ai = XiC + (d u u Q - x k u a>k - tu a )- h (dptfa - x k ip ajk - tip a )- 



du a ,i OLf^i 
DC 

+ (d^fap - Xk4>a/3,k ~ t(f) a p)— , (4.53) 

0(paP,i 

A 4 = tC + (d u u a - x k u a , k - tu a )^- + (d v (p a - x k tp ak - ttpa)-^- 

0U a dtp a 

dC 

+ (d<j>fap - Xk4>aP,k ~ t(p a p)—; • {4.54) 

OpaP 

In terms of force stress, pseudomoment stress and also with the canonical pseudo-momentum 
vector, the Poynting vector, the Eshelby stress tensor and the energy density the scaling 
flux densities become 

Yi := -Ai = xj Pji -tSi + d u o-ji uj + d 9 D j{ tpj - d^ H jH fa, (4.55) 
y :— A4 — XjVj — tH + d u pj uj + d^ D jt fa. (4.56) 

Eq. (I4.55P is the dilatational vector flux and (14.561) is some kind of generalized action. 
The balance law of scaling symmetry reads 



D t y - DiYi = (d u H — ) {piUi - u^a^) + (d v + -Jptfi - (d^ + 

+ (dtj, + —j—) {J ) >jO; i - ^ H ijk T ijk ) - (dp + ^-^jDijipij. (4.57) 



Eq. (I4.57P is the balance law for the scalar moment of momentum in the gauge theory 
of dislocations. It can be seen that we have two possibilities for the choice of d v and 
d<j). The first choice is close to micropolar and micromorphic elasticity (see, e.g., Lubarda 
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and Markenscoff [26j, Lazar [25|). Then, the scaling dimensions of the field variables are 



determined according 



d—2 . d d 

d u = — , dp = --, d^ = --, (4.58) 



where d = n+1 is the space-time dimension. In the present dynamic case for d = 4 = 3+1, 
we obtain d u — — 1, dp — — 2 and d^ — — 2. With this choice (14.581) . all the fields My, 
Uj, Mi and ifi have the same scaling dimension: — |. Thus, the minimal replacement 
( 12. 4p is scale invariant in the present case. Using (I4.58j) . the balance law (14.571) reads 

D t y-D i Y i = -2C di . (4.59) 

Thus, like in micromorphic elasticity the higher order terms break the scaling symmetry. 
In the present case and break the scaling symmetry. The corresponding M- 
integral is given by 

M (1) := / YimdS- [ DJW = 2 f C di dV. (4.60) 
Js Jv Jv 

The other choice for the scaling dimensions is 

d u = 2~ ' d v = ~~Y~' d<t> = ~~~2~' ^ 4 ' 61 ^ 

It is obvious that this choice is like in Maxwell's field theory of electromagnetic fields, 
where scalar and vector fields have the same dimensions, namely — (see, e.g., Felsager 
[27|). It can be seen that now the minimal replacement (12 .4[) is not scaling invariant 
because My has the dimension — | and 0^ has — ^=p. With d = 4, we obtain d u = — 1, 
= — 1 and dfy = —1. Now the scaling balance law (I4.57P reads 

D t y - = Pm - rr, ; 0; ; , (4.62) 

It is obvious that the source terms of the Euler-Lagrange equations (I2.16P and (12.171) 
break the scaling symmetry. The corresponding M-integral is now given by 

M (2) := / Y ini dS - I B t ydV = - I (pm - a.jOa) dV. (4.63) 
Js Jv Jv 



4.1.4 Gauge symmetry 

The gauge symmetry acts in the following way 

x'i = Xi, t' = t, u' a = u a - e f a , ip' a = ip a + e f a , (p' af3 = 4> a p + e f a ,p. (4.64) 
Then the components of the generator (13.111) take the form 

X ki = 0, T = 0, U a = -f a , ^ a = f a , $ a/3 = f a ,p. (4.65) 
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The density and flux current which can be derived from the infinite dimensional group 
generator H3.34f) by using the formulas (14. 9 p and (14.101) have the following form 

Gi := -Ai = -(Jjifj + Djifj - H jki f jtk , (4.66) 
Q' .1; - Pj fj + Djifo. (4.67) 

The corresponding conservation law reads 

D t g - DA = f jM H m = 0. (4.68) 

The global conservation law is 

/ Gi-riidS- [ D t gdV = 0. (4.69) 
Js Jv 

4.1.5 Addition of solutions 

The vector field v 6 is a generator of a divergence symmetry. The addition of solutions is 
given by 

x[ = Xi, t' = t, u' a = u a + eu a , (ff a = Lp a + e (p a , <p' af3 = <fr af3 + e 4> a p. (4.70) 

Then the components of the generator (13.111) take the form 

X ki = 0, r = 0, U a = u a , ^> a = (p a , $ Q/3 = Q/9 . (4.71) 

Using Betti's reciprocal theorem, the fields Bi and are of the form 

Bi = -u a a ai - ip a D ai + <p a p H a/3i , (4.72) 
B 4 = u a p a + 4> a p D af3 . (4.73) 

The notation <j ai , D ai and H a/ 3i means that u, ip and <f> are replaced by u, (p and <fi, 
respectively. Finally, the conserved fluxes are of the form 

A = -u a a ai - ip a D ai + 4> a/3 H a/3i + u a a ai + f a D ai - <p a/3 H a/3i , (4.74) 
A 4 = u a p a + (pap D a(3 -u a p a - 4> a/3 D a p. (4.75) 

The corresponding conservation law is a manifestation of Betti's reciprocal theorem for 



the gauge theory of dislocations and is a consequence of the linearity of A = [28|, |29 
In integral form the dynamical Betti reciprocal theorem is given by 

/ (u a &ai + <Pa Dai - <Paf3 H a f3i)rii dS - D t (u a Pa + <\>af3 D af3 ) dV = 

Js Jv 

/ (u a (?ai + Va D ai - (f) a/3 H a/3i )ni dS - D t (u a p a + <fi a/3 D aP ) dV. (4.76) 
Js Jv 
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4.2 Gauge- invariant currents 

Although all local conservation and balance laws for the total system can be brought 
in the gauge invariant form, this is not true for all the currents. Only the canonical 
currents (I4.18I) - (I4.21I) coming from the translational symmetry can be written in terms of 
gauge invariant quantities. For rewriting P ki , we use the Eq. (12.41) for the distortion, the 
dislocation current and torsion density tensor (12. 8p and the Eq. (12.221) . For V k we use also 
the same equations as for the canonical Eshelby-stress tensor, but instead of Eq. (12.221) . 
Eq. (12.211) is now needed. For rewriting Ti,, we use the Eq. (12 .4p for the kinematical 
velocity, the dislocation current density (12.81) and also the Eq. (12.211) . For Si we use the 
same equations as for the canonical energy density, but instead of (I2.2ip . Eq. (12.221) is 
used. The result reads as 

P kl = P® - D t (D ai <j) ak ) + D p (H api <f> ak ), (4.77) 

V k = VP -D p (D a p(j) ak ), (4.78) 

Si = Sf + D t (D m <p a ) - D^Hapi <p a ), (4.79) 

H = H®+D f ,(D aP <p a ), (4.80) 

where the gauge-invariant currents are defined by 

P k f := ~~ & $ki — &ai fiak + Dai lak ~ H aj3i T af3k , (4-81) 

V~l s) := -p a (3 ak + D aP T a(3k , (4.82) 
5j S := o a i v a — H aj 3i I aj 3, (4.83) 
ftte) ;= -C + Pa v a + D aP I aP . (4.84) 

From the canonical local conservation laws (14 . 1 8 j) — (14 . 2 1 j) we can obtain the gauge invariant 
conservation laws 

D t V { k s) - DiPg = 0, (4.85) 
D t H is) - DiS^ = 0, (4.86) 

since the divergence terms cancel each other out and 

D^tf^ (j) ak ) = 0, DiD^Hafx y a ) = 0, (4.87) 

because of the antisymmetry of H^j^ in the last two indices. By the help of the Gauss 
theorem, we get the momentum and energy conservation laws for the gauge invariant 
quantities in integral form 



4 g) : = / p lfm " / V t vi s) dV = 0, (4.88) 
Js Jv 

lis).- I s^riidS- [ D t H {g) dV = 0. (4.89) 
Js Jv 
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5 Configurational force, moment and power 



In this section we decompose the total Lagrangian into its elastic and dislocation parts. 
Both parts are coupled by the Euler-Lagrange equations. These coupling terms will appear 
as additional configurational forces, vector and scalar momenta. Consequently, the total 
system is conserved but they are not conserved separately. There is an exchange between 
the elastic and dislocation subsystems. 



5.1 Translations in space and time 

In the last section the conservation laws for the canonical pseudomomentum, material 
inertia vector and energy density have been calculated. These conservation laws exist 
only if the total Lagrangian C composed of the kinetic and potential energies of the 
elastic and the dislocation field part is taken into account. If we consider only one part 
of the field, the canonical quantities are not any more conserved and as a consequence 
configurational forces, moments and power appear. We give the currents for both parts of 
the total system, the elastic part C e and the dislocation part £di of the total Lagrangian 
density. 

5.1.1 Elastic part 

For the elastic part £ c we obtain from fj49l)- frCT0l) for the Eqs. fTCTSll - frOTTl the following 
canonical quantities of the elastic subsystem 

Pki := ~£c$ki ~ &ai u a,k, (5-1) 

V { k e) := -p a u a)k , (5.2) 

:= cr ai u a , (5.3) 

n {c) := -C e +p a ii a . (5.4) 

The local conservation laws (14.221) and (14.231) become now balance laws. At the right 
hand sides appear source terms called canonical configurational force and power. From 
Eqs. flSHD-flElD w e get 

D t vl c) -D t P^=F k , (5.5) 
D t H ic) - DiS^ = W, (5.6) 
where the configurational force and power densities read 

Fk = Pa <fia,k — Vaprf'aPMi (5-7) 
W = -p a <Pa + Cr al3 <j)al3- (5.8) 

The corresponding balance laws in integral form are given by 

4 6) := / PlfnidS- [ D t vi c) dV = - [ F k dV, (5.9) 
Js Jv Jv 

/( e ) := / S^mdS- [ D t H {c) dV=- [ WdV. (5.10) 
Js Jv Jv 
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A similar expression like ( 15. 7p was earlier given by Kadic and Edelen 12| . They called 
this force the 'elastic excess' force. 

For the elastic part C e of the total Lagrangian £ the gauge invariant currents (|4.8ip - 
( 14~84P are 



p( e -g) 
r ki 



s. 



(e-g) 
fc 

(e-g) 



■^-e &ki @ai Paki 
Pa ftaki 



&ai V a , 
^(e-g) ;= + 



Pa V 



From Eqs. (I5.1ip - (l5.14p . the balance laws can be written in the following form 



(e-g) 



w (g) , 



where 



7r(g) 
x k 

W (g) 



Pi lik @ij T-i 



ij -Lijkj 



@ij lij ; 



(5.11) 
(5.12) 
(5.13) 
(5.14) 

(5.15) 
(5.16) 

(5.17) 
(5.18) 



are the gauge invariant configurational force and power densities. In Eq. (15.171) we have 
two contributions to the configurational force density. The first term is the force density 
caused by a dislocation current density in presence of a momentum field pi (see also 
(30I . 31 . 32|). The second term is the Peach-Koehler force density [Hi] recovered in the 
gauge theory of dislocations. It gives the force density caused by a dislocation density 
in presence of a stress field a^. Here we call (I5.17P the dynamical Peach-Koehler force 
density. Eq. (15.181) is the power density lost from a moving dislocation to the medium 



33l . |30L |32J. The dynamical Peach-Koehler force is the analog of the Lorentz force of 
electromagnetism. If we multiply (I2.20p by Vi and the second equation of (12.91) by and 
add each other, we get 



ViPi + CTijPij - D^Vidi 



Cijlij. 



(5.19) 



Therefore, (15.181) is the rating density of the elastic material. The corresponding J and 
/-integrals read 



7(e-g) ._ 



j(e-g) ._ 



St s) nidS 



B t vt g) dV = - 



v 



H S) dV, 



D,H (c " g) dV = - W (s) dV. 



v 



(5.20) 
(5.21) 



v 



Static case 

For a static theory the Eshelby stress tensors of the elastic subsystem have the following 
form 



p(e) 
r H 

ki 



W e Ski — &ai ftak, 



(5.22) 
(5.23) 
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with the associated J-integrals 

( e ) _ / p( e )^ a Q — / -r( Sta ) 



Jr = / P£ } nidS= / H / a^dV, (5.24) 

Js Jv Jv 

4 e ' g) = I P^ g) nidS= [ ^ K) dV = - [ a ij T ijk dV = [ e kjl a ijau dV, (5.25) 
Js Jv Jv Jv 

where Tij k = eijkOiu with the usual dislocation density tensor an. The gauge invariant 
Eshelby stress tensor (15.231) is in agreement with the Eshelby stress tensor in incompatible 
elasticity derived by Kirchner [341 ]. Lazar and Kirchner 35j] . It can be seen that the 
configurational force density in (15.241) is just given in terms of the force stress tensor and 
the gradient of the plastic distortion. But this is not the correct Peach-Koehler force. On 
the other hand, the gauge invariant expression (15.251) gives the correct gauge-invariant 
Peach-Koehler force. 

5.1.2 Dislocation part 

In the same manner we obtain for the canonical currents of the dislocation part 

Phi '■= —£di Ski — D ai ip a)k + H aj3i (j) a p,ki (5.26) 

Pf :=-D af}( f> a p, k , (5.27) 

Si '■= D ai tp a - H al3i <f> af3 , (5.28) 

H {A) := -C di + D a p<p af3 (5.29) 

with the J and /-integrals of the dislocation subsystem 

4 d) := / P^n.dS- [ D t V ( k d) dV= I F k dV, (5.30) 
Js Jv Jv 

j( d ) : = [ sfrndS- [ D t n (d) dV= [ WdV. (5.31) 
Js Jv Jv 

The gauge invariant currents of the dislocation part C^i of the total Lagrangian £ read 

P k i ^ '■= — Aii ^ki + D ai I ak — H a pi T a p k , (5.32) 

'P < k^''- = P ) apT a i 3k , (5.33) 
S\ d ^ '■= —H a /3il a f3, (5.34) 
^ :=-£ di + D a/ ,V (5.35) 

From them the configurational force and power can be calculated in the same way, as 
we did for the elastic part. If we use the Eqs. (I5.26I) - (I5.29I) and (I5.32I) - (I5.35I) for the 
calculation of the canonical and gauge invariant force and power, the same results as in 
the Eqs. f )5.15p . (15.161) and (15.171) . (I5.18P appear, but with opposite signs. The result is 

Jt e) ■= [ P£ s) nidS- [ D t vt e) dV= [ J#>dV, (5.36) 
Js Jv Jv 

j(d-g) ;= I s^mdS - [ D t H (d - g) dV= [ W is) dV. (5.37) 
Js Jv Jv 
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If we multiply (12. 22ft by 1^ and the second equation of (12.101) by Hij k and add each other, 
we obtain 

lijDij + HijkTijk + Dk(IijHijk) — (5.38) 

This is the power density lost by moving dislocations. 

Finally, the J and the I integrals of the total system may be decomposed into its 
elastic and dislocation parts according to 

J k = 4 e) + jf = 0, (5.39) 
/ = /( c ) + /( d ) = 0. (5.40) 

The configurational force and power densities appearing in the elastic and dislocation 
subsystems are the result of the interaction of both subsystems. Without interaction 
between the elastic and dislocation parts no configurational force and power densities 
would appear and J^ e ' = 0, = J^ d ' = 0, 1^ = 0. Of course the same is true for the 
gauge invariant parts. 

Static case 

The static Eshelby stress tensors of the dislocation subsystem are 

P kf = w di Ski + H a/ 3i 4>af3,k = W di 5 ki - e m H a i4> a ^ k , (5.41) 

Pki ^ = Wdi ~ Haf3i T a /3k = Wdi $ki + ^if3lH a iT a p k , (5-42) 



with the J integrals 



4 d) = / P$nidS=- / ^ sta) dy = - / VijfakdV, (5.43) 
Js Jv Jv 

Jt &) = [ P£ S) ni dS = -f 4 PK) dV= [ a t3 T t3k dV. (5.44) 
Js Jv Jv 

For a dislocation theory without force stresses, cr^ = 0, the Eshelby stress tensors 
(I5.4ip and f!5.42j) are divergenceless and the J integrals (I5.43P and (15.441) are zero. This 
is the case for a (force) stress-free configuration with pseudo-moment stresses only. Then 
the Euler-Lagrange equation (I2.22p takes the form 

DiHafu = 0. (5.45) 



If we use the so-called Einstein choice C2 = —c\ and C3 = — 2ci [36|, [37|, Eq. (15.451) (or 



Eq. (I2.30p ) reduces to the incompatibility condition for the elastic strain 22] : ince = 



where e = |(/3+/3 T ). Anyway, the Eshelby stress tensors (15.411) and (I5.42p are the Eshelby 



stress tensors of the dislocation part. Recently, Li et al. 38j and Li 39j have used such an 
Eshelby stress tensor of dislocations without interaction with a force stress tensor. They 
called it compatibility momentum tensor. Furthermore, they have found conservation 
laws associated with compatibility conditions of continua. But, in fact, they investigated 
the static dislocation part without force stresses and they found ( 15.451) as compatibility 
condition and the J^ d \ and integrals of the dislocation subsystem. It is obvious 
that such conservation laws are not a "new" class of conservation laws. If the Peach- 
Koehler force is zero, they are conservation laws. Thus, it is misleading to call (15.411) a 
compatibility momentum tensor. 
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5.2 Rotations in space 

5.2.1 Elastic part 

In the same way as we did for the translations, we find for the canonical currents for the 
elastic part of the system described by £ e i 

M ti ■ = tkrnj (x m P\f + u m a 3i ) , (5.46) 
Mf : = e kmj (x m Vf + u m Pj ), (5.47) 

where Pjjf and P- are given by (15.11) and (15. 2p . respectively. The rotational balance law 
of the elastic part is calculated 

D t M$ - DiM^f = 6 kmj (x m J 7 j + (v m Pj - /3 ml aji - /3j m <7y) - Lp m Pj + (j) m l^jl + famuli), 

(5.48) 

where so-called configurational vector moments appear as source terms. In Eq. (15.481) 
the first part is a moment produced by the configurational force J-j, the next three parts 
are the isotropy condition which vanishes if the constitutive relation between cr and f3 is 
isotropic. The other source terms are caused by the gauge potential which are the plastic 
fields. The L integral reads now as 

L k c) := / M^friidS - / D t M k c) dV = / e k j m {x m Tj + (v m pj - f3 m i°ji ~ Pim^ij) 
Js Jv Jv 

- PmPj + <f>ml(7jl + fylwPlj) & V - ( 5 - 49 ) 

Static case 

In the static case, the angular momentum tensor (15.461) is given in terms of the static 
Eshelby stress tensor (I5.22p . Eventually, the static L-integral of the elastic subsystem 
reads 

L k e) = / M^frii dS = e k j m (x m ^f t&) - {fi m i<Jji + PimVij) + (fimiVji + 4>i m crij) dV. 
Js Jv 

(5.50) 

5.2.2 Dislocation part 

For the dislocation part of the system with the Lagrangian and Pjf , given 
by (I5.26P and (15.271) we get for the canonical currents 

Mj^f : = e km j (x m Pji + Dji <Pm - Hjii <Pmi - H t ji(j)i m ), (5.51) 
M ( k d) : = e kmj (x m Vf ] + D fl 4> ml + Dij Im ). (5.52) 

Now the rotational balance law of the dislocation subsystem is given by 

D t M£ } - D^M^f = -e fcmj (x m T$ - (v m Pj - f3 m i(?ji - P>i m o Xj ) + u m pj + u m ,iOji + u Iim oy 

— ImlDjl — IlmDlj + - TmilHjU + Tu m Huj). (5.53) 
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The corresponding L integral reads 

= J M^frii dS - J T) t Mf ] dV = - J ekjm^XmFj ~ (VmPj - PmlVjl ~ Plm&lj) 

+ u mPj + u mm + Ul , m a, - I mlDjl - IimDij + i T milHjil + T lim H Uj ) dV. (5.54) 

Finally, we get the decomposition of the L integral into the elastic and dislocation parts 
according 

L k = Lt ] + Lf\ (5.55) 

Static case 

For the static case, the angular momentum tensor of the dislocation part is 

= € kmj {XmPji^ — HjU(f) m i — Hiji(pi m ), 

= tkmj { x m,Pji ^ — QinHj n <p m l) — 5 ki Hi n (f)i n + (5.56) 

with the static L integrals as follows 

Lf ) = J Mlfrii dS = J t kjm (^X m jj Sta) - (PmlVjl + Plrn^lj) + «m,/^7 + Ul,m^lj 

+ - T mil H jU + T Um H u ^j dV. (5.57) 

Of course, for an isotropic material in a force stress-free state Eq. (I5.57P is a conservation 
law. 

5.3 Scaling in space 
5.3.1 Elastic part 

For the scaling we can also calculate for the Lagrangian C e the canonical currents 

y/ e) := Xj - t + d u crji Uj , (5.58) 
3#) ■= Xj vf - t H {c) + d uPj Uj , (5.59) 



where P$\ vf \ Sf and are given by fl5.ip - fl5.4l) . The corresponding balance law 
reads 

D t y^ - DtY^ = -tW + x k F k + ^ Pa( p a - 1 (Tapfa. (5.60) 
For the path-independent M k integral we obtain 

M®:= j Y^mdS- j D t y^dV = J (tW - x h T k -^p a (j> ol + ^a ol p ( t> ot p)dV. 

(5.61) 

22 



Static case 

In the static case, the scaling vector (15.581) is given in terms of the static Eshelby stress 
tensor f 1 5 . 2 2 j) as follows 

y/ e) = Xj Pj° + d u crji Uj . (5.62) 
The static M integral of the elastic subsystem reads 

^ = Is Y ' e)nidS = ~Jv fa*" ~t aa ^) dV - (5 ' 63) 
5.3.2 Dislocation part 

For the canonical currents of the dislocation part from Eqs. (I5.26f) — ( T5.29H we get 

y/ d) := Xj - t Sf } + d v Dji i Pj - d^fa, (5.64) 
y (d) := Xj Vf - t H W + Dji fa. (5.65) 

In general, the balance law derived from the dislocation part is given by 

D t y {d) - DiF/ d) = tW - x k F k + d^paipa - d^a^ap - (d v + ^—^)D aj3 ip a3 

+ (dj, + ^-y^) A*/30 Q /3 ~\{ d <t> + —Y~) Ha f 3kTa f 3k - (5.66) 
Using (14.581) . the balance law (I5.66P gets the form 

D^ (d) - D,y/ d) =tW-x k F k - ^ Pa <p a + ~a af 3<f) aP - 2C di . (5.67) 
In integral form it reads 

:= J Y^mdS- j D t y^dV = -J (tW -x k F k -^p a( j) a + ^a a p(j )al3 -2jC, d ?)dV. 

(5.68) 

If we use (I4.6ip . we obtain from (15.661) the following balance law 

D^ d ) - DiY} d) = tW-x k F k - d -^- Pa <p a + d -^- a af3 <p af3 . (5.69) 

For vanishing source terms of the Euler-Lagrange equations (12.211) and (l2.22l) . the balance 
law (15.691) reduces to a conservation law. In integral form we have 

M<™ := J Y^mdS- J D t y^dV = - J (tW-x k F k -^p a ct> a + ^a aP <f> a pyV. 

(5.70) 

Finally, the decomposition of the M integral into the elastic and dislocation parts read 

M (l) = M (e) +M (l-d) ; ( 571 ) 
M (2) = M (e) +M (2-d)_ 
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Static case 

In the static case, the scaling vector (15.641) is given in terms of the static Eshelby stress 
tensor (I5.4ip as follows 

Y^ d = xj Pjf + Hjn (f)ji = Xj — d,/, eu m Hj m </> 3 -j. (5.73) 

The static M integrals of the dislocation subsystem with the choices (14.611) and (14.611) are 
given by 

M (i-d) = J F .(d) n . dS = J ^ ^(sta) _ n + dy . (5 

= j Y^n, dS = J (x k a) - ^ a^p) dV. (5.75) 
If cry = (stress- free state), then only ( 15.75ft is a conservation law. 

5.4 Gauge symmetry 

5.4.1 Elastic part 

With the following currents 

Gf ■= -Vjifj, (5.76) 

G ie) := -Pjfj, (5.77) 

we get 

D t g® - D.Gf = - Pj f d + ajifo. (5.78) 

It is obvious that for fi = constant, Eq. (15.781) becomes a conservation law. 
For the static case, we obtain 

DiGV = -ajifai. (5.79) 

5.4.2 Dislocation part 

For the dislocation part, the flux quantities are 

:= Djifj — Hjkifj^, (5.80) 

:= Drf^ (5.81) 

with the balance law 

D t g^ - DiGf = p 3 f 3 - o,,f,,. (5.82) 

If fi = constant, Eq. (15.821) becomes a conservation law. 
For the static case, we get 

G\ ^ = —Hjkifj,k = £ikiHjifj t k, (5.83) 



24 



and the balance law 

.(d) 



D i Gt ) = ajif^. (5.84) 

For Uji = (stress- free state), the balance law (15.841) becomes a conservation law. This is 
the case for a pure dislocation theory without force stresses. If fj^ = constant, we recover 



the symmetry of constant pre-distortion used by Li et al. [38| and Li 39 



5.5 Addition of solutions 

5.5.1 Elastic part 

From Eqs. (I4.74p and (I4.75P we obtain for the elastic subsystem 

A (C) = -u a v a i + u a a ai , (5.85) 
Af> =u a p a - u a p a . (5.86) 

They fulfill the corresponding balance law 

BtA^ + DiAf ] =p a ti a - p a U a - Oat U a ,i + Oai u a ,i (5.87) 

For the static case we have the balance law 

DiAf' = -Oai U a ,i + Oai U a ,i- (5.88) 

5.5.2 Dislocation part 

From Eqs. (I4.74p and (I4.75P we obtain for the dislocation subsystem 

Af> = -<f a D ai + <Paf3 H api + Lp a Am ~ 4>aP H a (3i, (5.89) 

4 d) =4>apD aP -<l> a pD a p. (5.90) 
These quantities lead to the following balance law 

D t Af } + DiAf =p a y a - pa V?o - &ai 4>m + °ai 4>ai- (5.91) 

If the source terms are zero (p a = 0, p a = 0, a a % = and a ai = 0) the balance law (15.911) 
becomes a conservation law. 
In the static case we have 

= 4>af3 H a! 3i — 4>al3 H a j3i = Cifik (0a/3 H ak — <j) a p H a k) (5.92) 

with the balance law 

DiAf } = -Oai 4>ai + Oai <fiai- (5.93) 

If Oai = and a Q i = (stress-free state), the balance law (I5.93P becomes a conservation 
law. 



25 



6 Conclusion 

We have presented the equations of motion for the translational gauge field theory of 
dislocations with asymmetric stresses. For this purpose we have chosen the most general 
linear isotropic constitutive equations. The Lie-point symmetries leaving these system of 
equations form-invariant have been discussed. According to the Noether theorem, from 
the variational and divergence symmetries for the total Lagrangian density the currents 
in canonical form and the local conservation laws for the translational and rotational 
invariance have been derived. The form of the local balance law for the broken scaling 
symmetry has been found. The gauge-invariant currents have also been obtained. Taking 
into account only the elastic or dislocation part of the total Lagrangian density we have 
shown how the local conservation laws for the continuous transformation of translation and 
rotation turned into balance laws giving rise to a configurational force and moment. In this 
manner we have found the dynamical Peach-Koehler force. In addition we have calculated 
the configurational power for both parts. Since no external forces and moments act on the 
whole system described by the total Lagrangian the energy, linear and angular momentum 
are conserved. We also derived the conservation and balance laws for the gauge symmetry 
and the addition of solutions. Using the divergence symmetry of addition of solutions, we 
were able to derive a reciprocity theorem for the gauge theory of dislocations. In addition, 
we have derived the conservation laws for (force) stress-free states of dislocations . If we 
identify 0^ = — with the plastic distortion, our results and the formal structure of 
currents keep valid in gradient plasticity. 
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